pH Regulates Ligand Binding to an Enzyme Active Site by Modulating Intermediate Populations.

Understanding the mechanism of ligands binding to their protein targets and the influence of various factors governing the binding thermodynamics is essential for rational drug design. The solution pH is one of the critical factors that can influence ligand binding to a protein cavity, especially in enzymes whose function is sensitive to the pH. Using computer simulations, we studied the pH effect on the binding of a guanidinium ion (Gdm+) to the active site of hen egg-white lysozyme (HEWL). HEWL serves as a model system for enzymes with two acidic residues in the active site and ligands with Gdm+ moieties, which can bind to the active sites of such enzymes and are present in several approved drugs treating various disorders. The computed free energy surface (FES) shows that Gdm+ binds to the HEWL active site using two dominant binding pathways populating multiple intermediates. We show that the residues close to the active site that can anchor the ligand could play a critical role in ligand binding. Using a Markov state model, we quantified the lifetimes and kinetic pathways connecting the different states in the FES. The protonation and deprotonation of the acidic residues in the active site in response to the pH change strongly influence the Gdm+ binding. There is a sharp jump in the ligand-binding rate constant when the pH approaches the largest pKa of the acidic residue present in the active site. The simulations reveal that, at most, three Gdm+ can bind at the active site, with the Gdm+ bound in the cavity of the active site acting as a scaffold for the other two Gdm+ ions binding. These results can aid in providing greater insights into designing novel molecules containing Gdm+ moieties that can have high binding affinities to inhibit the function of enzymes with acidic residues in their active site.

Cosolvent effects on the growth of amyloid fibrils.

Cells are equipped with cosolvents that modulate protein folding and aggregation to withstand water stress. The effect of cosolvents on the aggregation rates depends on whether the polypeptide sequence is an intrinsically disordered protein (IDP) or can fold into a specific native structure. Cosolvents, which act as denaturants generally slow down aggregation in IDPs, while expediting it in globular proteins. In contrast, protecting osmolytes facilitate aggregation in IDPs, while slowing it down in globular proteins. In this review we highlight the recent computational approaches to gain insight into the role of cosolvents on the aggregation mechanism of IDPs and globular proteins. Computer simulations using the molecular transfer model, which implements the cosolvent effects in coarse-grained protein models in conjunction with enhanced sampling techniques played an important role in elucidating the effect of cosolvents on the growth of amyloid fibrils.

Role of Guanidinium-Carboxylate Ion Interaction in Enzyme Inhibition with Implications for Drug Design.

Guanidinium cation (Gdm+) interacts strongly with amino acids of different polarities modulating protein structure and function. Using density functional theory calculations and molecular dynamics simulations, we studied the interaction of Gdm+ with carboxylate ions mimicking its interaction with acidic amino acids and explored its effect in enzymatic folding and activity. We show that, in low concentrations, Gdm+ stabilizes carboxylate ion dimers by acting as a bridge between them, thereby reducing the electrostatic repulsion. We further show that this carboxylate-Gdm+-carboxylate interaction can have an effect on the structure-activity relationship in enzymes with active sites containing two acidic residues. Using five enzymes (hen egg white lysozyme, T4 lysozyme, HIV-1 protease, pepsin, and creatine kinase), which have two acidic amino acids in their active sites, we show that, in low concentrations (<0.5 M), Gdm+ strongly binds to the enzyme active site, thereby potentially inhibiting its activity without unfolding it. This can lead to misleading conclusions in experiments, which infer the extent of enzyme unfolding from activity measurements. However, the carboxylate-Gdm+-carboxylate specific interaction can be exploited in drug discovery as drugs based on guanidinium derivatives are already being used to treat various maladies related to muscle weakness, cancer, diabetes etc. Guanidinium derivatives can be designed as potential drug molecules to inhibit activity or functioning of enzymes, which have binding pockets with two acidic residues in close vicinity.

Role of Disulfide Bonds and Topological Frustration in the Kinetic Partitioning of Lysozyme Folding Pathways.

Disulfide bonds in proteins can strongly influence the folding pathways by constraining the conformational space. Lysozyme has four disulfide bonds and is widely studied for its antibacterial properties. Experiments on lysozyme infer that the protein folds through a fast and a slow pathway. However, the reasons for the kinetic partitioning in the folding pathways are not completely clear. Using a coarse-grained protein model and simulations, we show that two out of the four disulfide bonds, which are present in the α-domain of lysozyme, are responsible for the slow folding pathway. In this pathway, a kinetically trapped intermediate state, which is close to the native state, is populated. In this state, the orientations of α-helices present in the α-domain are misaligned relative to each other. The protein in this state has to partially unfold by breaking down the interhelical contacts between the misaligned helices to fold to the native state. However, the topological constraints due to the two disulfide bonds present in the α-domain make the protein less flexible, and it is trapped in this conformation for hundreds of milliseconds. On disabling these disulfide bonds, we find that the kinetically trapped intermediate state and the slow folding pathway disappear. Simulations mimicking the folding of protein without disulfide bonds under oxidative conditions show that the native disulfide bonds are formed as the protein folds, indicating that folding guides the formation of disulfide bonds. The sequence of formation of the disulfide bonds is Cys64-Cys80 → Cys76-Cys94 → Cys30-Cys115 → Cys6-Cys127. Any disulfide bond that forms before its precursor in the sequence has to break and follow the sequence for the protein to fold. These results show that lysozyme also serves as a very good model system to probe the role of disulfide bonds and topological frustration in protein folding. The predictions from the simulations can be verified by single-molecule fluorescence resonance energy transfer or single-molecule pulling experiments, which can probe heterogeneity in the folding pathways.

Salt Effects on Protein Folding Thermodynamics.

Salts differ in their ability to stabilize protein conformations, thereby affecting the thermodynamics and kinetics of protein folding. We developed a coarse-grained protein model that can predict salt-induced changes in protein properties by using the transfer free-energy data of various chemical groups from water to salt solutions. Using this model and molecular dynamics simulations, we probed the effect of seven different salts on the folding thermodynamics of the DNA binding domain of lac repressor protein ( lac-DBD) and N-terminal domain of ribosomal protein (NTL9). We show that a salt can act as a protein stabilizing or destabilizing agent depending on the protein sequence and folded state topology. The computed thermodynamic properties, especially the m values for various salts, which reveal the relative ability of a salt to stabilize the protein folded state, are in quantitative agreement with the experimentally measured values. The computations show that the degree of protein compaction in the denatured ensemble strongly depends on the salt identity, and for the same variation in salt concentration, the compaction in the protein dimensions varies from ∼4% to ∼30% depending on the salt. The transition-state ensemble (TSE) of lac-DBD is homogeneous and polarized, while the TSE of NTL9 is heterogeneous and diffusive. Salts induce subtle structural changes in the TSE that are in agreement with Hammond's postulate. The barrier to protein folding tends to disappear in the presence of moderate concentrations (∼3-4 m) of strongly stabilizing salts.

Contrasting Effects of Guanidinium Chloride and Urea on the Activity and Unfolding of Lysozyme.

Cosolvents play an important role in regulating the stability and function of proteins present in the cell. We studied the role of cosolvents, urea and guanidinium chloride (GdmCl), which act as protein denaturants, in the catalytic activity and structural stability of the protein lysozyme using activity measurements, spectroscopy, and molecular dynamics simulations. We find that the activity of lysozyme increases on the addition of urea, whereas it decreases sharply on the addition of GdmCl. At low GdmCl concentrations ([GdmCl] < 4 M), the activity of lysozyme decreases, even though there is no significant perturbation in the structure of the lysozyme folded state. We find that this is due to the strong interaction of the Gdm+ ion with the residues Asp52 and Glu35, which are present in the lysozyme catalytic site. In contrast, urea interacts with Trp63 present in the loop region present near the active site of lysozyme, inducing minor conformational changes in lysozyme, which can increase the activity of lysozyme. At higher denaturant concentrations, experiments show that GdmCl completely denatures the protein, whereas the folded state is stable in the presence of urea. We further show that GdmCl denatures lysozyme with the disulfide bonds intact in the protein, whereas urea denatures the protein only when the disulfide bonds are broken using reducing agents.

Osmolyte Effects on the Growth of Amyloid Fibrils.

Understanding the role of naturally occurring protective osmolytes, such as trimethylamine N-oxide (TMAO), in the growth of amyloid fibrils implicated in neurodegenerative diseases is important to prevent fibril growth. The effect of TMAO on the growth of amyloid fibrils formed by the Sup35 prion peptide NNQQNY is studied using molecular dynamics simulations. The free-energy surface for the growth of the protofibril shows three major basins, corresponding to the free state where the peptide is in solution, the docked state where the peptide in solution interacts with the surface of the protofibril, and the locked state where the peptide is tightly bound to the protofibril, becoming a part of the fibril. The free-energy surface in the presence of TMAO shows that TMAO stabilizes the locked state of the peptide compared with the free state, indicating that TMAO aids in fibril growth. Dissecting the interaction of TMAO with individual amino acids in the peptide shows that TMAO interacts both directly and indirectly with the amino acids, depending on the nature of the side chains. The methyl groups in TMAO interact strongly with the hydrophobic aromatic ring in the Tyr residue. In the locked state, the surface area of Tyr available for interaction with TMAO decreases; as a result, the Tyr residue in the peptide flips out from the locked position, increasing the fluctuations of the peptide locked in the protofibril. Such strong direct interactions of amino acids with TMAO destabilize the folded or aggregated states of proteins. The overall increased stability of the peptide locked in the protofibril by TMAO is due to entropic or indirect interactions with the backbone Asn and Gln residues, which form major components of the NNQQNY peptide.